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We present Fabry-Perot absorption-line spectroscopy of the SBO galaxy NGC 
7079. This is the first use of Fabry-Perot techniques to measure the two- 



dimensional stellar kinematics of an early-type disk galaxy. We scan the infrared 
Call line using the Rutgers Fabry-Perot (RFP), to obtain kinematic data extend- 
ing to /-band surface brightness jij ~ 21 mag arcsec -2 , in a field of radius ~ 40". 



The kinematic data, consisting of line-of-sight velocities and velocity dispersions, 
are in good agreement with data obtained along the major axis of the disk with 
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standard slit spectroscopy. Comparison of the exposure times required for slit 
and RFP spectroscopy to reach the same limiting magnitude shows that the RFP 
is significantly more efficient for mapping absorption-line galaxy kinematics. We 



use the velocity data, together with our own deep broad-band photometry, to 
measure the bar pattern speed, Q, p , of NGC 7079 with the model- independent 
Tremaine- Weinberg (TW) method. We find Q p = 8.4 ± 0.2 km s _1 arcsec _1 ; this 
is the best-constrained pattern speed ever measured for a bar using the TW 
method. From the rotation curve, corrected for asymmetric drift, we calculate 
the co-rotation radius and find that the bar ends just inside this radius. The 
two-dimensional character of these data allow us to show that the TW method 
is sensitive to errors in the position angle (PA) of the disk. For example, a PA 
error of 2° can give errors ~ ±25% in Q p . 
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1. Introduction 

Recent experience has demonstrated the utility of integral-field spectroscopy to studies 
of galaxy dynamics. Verolme et al. (2002), comparing mass models of M32 constrained by 
integral-field spectroscopy from SAURON (Bacon et al. 2001) and models constrained only 
by data along 4 slits extracted from the full two-dimensional (2-D) data set, found that the 
mass model parameters, including black hole mass, mass-to-light ratio and, in particular, 
inclination were substantially better constrained with the full 2-D data than with just the 
slits. 

One approach to integral-field spectroscopy is the use of an imaging Fabry-Perot inter- 
ferometer (FP), essentially a tunable, narrow-band filter. The Rutgers Fabry-Perot (RFP), 
formerly 4 at CTIO 5 , found considerable application to the emission-line study of kinematics 
and mass modeling of disk galaxies {e.g. Buta & Purcell 1998; Buta et al. 1999 & 2001; 
Beauvais & Bothun 1999 & 2001; Palunas & Williams 2000; Weiner et al. 2001 a,b). How- 
ever, FP stellar absorption-line spectroscopy in galaxies is significantly more difficult due to 
the lower surface brightness levels, the shallowness and width of the spectral features, and 
the effects of time-variable observing conditions on reconstructing line profiles necessarily 
sampled sequentially over an extended period of observation. Because of these difficulties, 
FP stellar absorption-line spectroscopy has not been previously successfully attempted in 
galaxies, although absorption-line spectroscopy of individual stars in globular clusters has 
been accomplished with considerable success (Gebhardt et al. 1994, 1995 & 1997). In this 
paper, we study the kinematics of a bright, barred, disk galaxy with stellar absorption-line 



1 Based in part on observations carried out at the European Southern Observatory (Prop. No. B-0329) 

4 The RFP has since been de-commissioned. Replacement by the next generation RFP is scheduled for 
2004. 

5 The Cerro Tololo Inter- American Observatory (CTIO) is part of The National Optical Astronomy Ob- 
servatories (NOAO) which are operated by the Association of Universities for Research in Astronomy (AURA 
Inc.), under cooperative agreement with the National Science Foundation. 
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spectroscopy using the RFP, which is well-suited to galactic spectroscopy because of its large 
field-of-view, high throughput, and appropriate spectral resolution (R = \/5\ ~ 2000). 

More than half of all high surface brightness disk galaxies are barred (Knapen 1999; 
Eskridge et al. 2000). The fundamental parameter in the dynamics of barred galaxies is the 
pattern speed of the bar, Q p . This is usually expressed in terms of the quantity TZ = Di/as', 
here as is the length of the semi- major axis of the bar and Dl is the Lagrangian/corotation 
radius, where, in the bar's rest-frame, the centrifugal and gravitational forces balance. A 
self-consistent bar with orbits aligned with the bar has 1Z > 1 (Contopoulos 1980); when 
1.0 < 1Z ^ 1.4, a bar is termed fast. 

Because bars have strong quadrupole moments, they lose angular momentum efficiently 
in the presence of a dense dark matter (DM) halo, as first shown by perturbation theory cal- 
culations (Weinberg 1985). Early simplified simulations (Little k Carlberg 1991, Hernquist 
k Weinberg 1992) found agreement with Weinberg's calculations. The presence of DM in 
disk galaxies is required by their rotation curves, which stay flat out to large radii (Rubin et 
al. 1980; Bosma 1981). Determining the relative amounts of luminous and dark matter in 
the region of the optical disk, however, is a difficult problem, due to degeneracy in the con- 
tributions of the two components to the rotation of axi-symmetric systems. The maximum 
disk hypothesis (van Albada k Sancisi 1986), which requires the maximum amount of disk 
matter consistent with the inner rotation curve, often accounts for the shapes of the rotation 
curves of high surface brightness (HSB) galaxies (e.g. Kalnajs 1983; Kent 1986; Corsini et 
al. 1999; Palunas k Williams 2000). Others studies, however, have argued that maximum 
disks cannot account for other observational properties (e.g. van der Kruit 1995; Bottema 
1997; Courteau k Rix 1999, Kranz et al. 2001). Debattista k Sellwood (1998) realized, from 
fully self-consistent simulations, that fast bars require that barred galaxies have maximum 
discs. At present, there is a lively debate about this suggestion (Tremaine k Ostriker 1999; 
Debattista k Sellwood 2000; Weinberg k Katz 2002; Valenzuela k Klypin 2003; O'Neill k 
Dubinski 2003). 

One of the issues surrounding this debate concerns just how rapidly bars rotate. The 
number of barred galaxies in which 1Z has been measured is still relatively small. Evidence 
that bars are fast comes mostly from simulations of the hydrodynamics of gas in the bar 
region (e.g. Athanassoula 1992; Lindblad et al. 1996). The example of NGC 4123 (Weiner et 
al. 2001b) is especially instructive because the hydrodynamical simulations require both a fast 
bar and a maximum disk to match the observed gas shocks. A direct (model-independent) 
method for measuring bar pattern speeds was devised by Tremaine k Weinberg (1984). The 
Tremaine- Weinberg (TW) method is contained in the following simple equation 



XQ v sini = V. 



(1) 
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Here, X = J h(Y) X Z dX dY , V = J h(Y) V los ZdXdY, £ is the surface brightness, h(Y) 
is an arbitrary weighting function, Vi os is the line-of-sight velocity, % is the disk inclination 
(throughout, we use the convention % — for face-on) and (X, Y) are galaxy-centered coordi- 
nates measured along the disc's apparent major and minor axes, respectively. Although the 
integrals in X and V range over — oo < X < oo, if the disk is axisymmetric at \X\ > X a , then 
integration over — X Q < X < X Q suffices. The integration in Y is over an arbitrary range, 
because of h(Y); for example a weight function proportional to a delta function corresponds 
to a slit parallel to the disk major axis. For a number of such slits, plotting V against X gives 
a straight line with slope fl p sini The TW method has previously been applied to a handful 
of galaxies, all with slit data (Merrifield & Kuijken 1995; Gerssen et al. 1999; Debattista et 
al. 2002; Aguerri et al. 2003; Gerssen 2003). In all cases, bars were found to be consistent 
with 1.0 < 1Z ^ 1.4. However, TW measurements are technically challenging, and in some 
cases large uncertainties result. Thus, for example, Valenzuela & Klypin (2003) argue that 
R = 1.7, which they find in their simulations, is still consistent with observations. While 
correct, this claim overlooks the fact that this limit is reached only in those galaxies with 
the largest uncertainties (Aguerri et al. 2003, hereafter ADC03) and ignores other sources of 
scatter (Debattista 2003, hereafter D03). There exists, therefore, a need for at least a few 
high precision measurements of 71. 

Integral field velocity data permit TW measurements of £l p with greater precision, be- 
cause of the improved spatial coverage. In this paper, we present observations of NGC 7079 
using the RFP interferometer, which allowed us to measure 2-D kinematics down to a surface 
brightness /!/ ~ 21 mag arcsec -2 . In §2 we describe NGC 7079, emphasizing those properties 
which make it a suitable system for this analysis. In §3, we present photometry of NGC 7079. 
This is the first use of the RFP for absorption-line spectroscopy of galaxies, and therefore 
in §4 we present the spectroscopy and reduction methods. We show that the resulting data 
are in good agreement with data obtained by standard slit spectroscopy. In §5, we study 
the kinematics of NGC 7079 and measure Q p with the TW method. We show at this point 
that the measured f2 p is quite sensitive to small errors in the position angle (PA) of the disk. 
Finally, in §6, we discuss our results and conclusions. 
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2. NGC 7079 

In choosing candidates for observation, SBO galaxies were selected from various cata- 
logues and examined in the Digitized Sky Survey 6 . Candidates were chosen with a bar at 
a position angle intermediate between the minor and major axes of the disk (as required 
by the TW method), a number of bright stars within 80" of the galaxy center (to permit 
determination of temporal atmospheric transmission variations between the RFP frames) 
and lacking strong spirals, patchy dust and significant companions. The galaxies needed to 
be bright and large, and have known recession velocities (required for efficient use of the 
RFP). Finally, we chose to use the Call 8542.09 A absorption line because it is quite strong, 
and the combination of the red color of SBO galaxies and the moderately high CCD detector 
quantum efficiency in this portion of the spectrum maximizes the measurement efficiency 
here. Since this region of the spectrum contains a host of night sky emission features (see 
Osterbrock et al. 1996), we selected candidate galaxies with red-shifts that placed the Call 
line in relatively emission-free regions. 

NGC 7079 (ESO 287-36) satisfied all these conditions. NGC 7079 is classified as 
(L)SB(r)0° (Buta 1995), is large and bright (2f3 x B T = 12.46 [de Vaucouleurs et 
al. 1991, hereafter RC3]) and has a measured recession velocity (2670 ±25 km s _1 , [Lauberts 
& Valentijn 1989]). The red-shifted wavelength of the Call line is therefore 8618A; assum- 
ing a rotation velocity of about 260 km s _1 (Bettoni & Galletta 1997) and allowing for the 
width of the Call stellar line and for the gradient in wavelength across the RFP frame, we 
needed to scan the galaxy from 8608A to 8631 A. The sky has no emission between 8598 A 
and 8613A and some weak O2 emission (which can be subtracted out) between 8613A and 
8631 A. Bettoni & Galletta (1997) detected 0[III] A = 5007A emission extending out to 
~ 15", but this gas is faintly emitting and is not expected to affect measurement of Q p . 
NGC 7079 is the brightest member of a group containing seven galaxies (Garcia 1993) but 
the nearest large neighbor, ESO 287-37, is over 0?5 away from it, or 363 kpc (in projection) 
at a distance of 38 Mpc (we assume throughout this paper that H = 70 km s _1 Mpc _1 ). 
The other group members are more than twice as far away, so we expect that NGC 7079 
should not be significantly perturbed. 



6 The Digitized Sky Survey was produced at the Space Telescope Science Institute (STScI) and funded by 
the National Aeronautics and Space Administration. 
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3. Surface Photometry 

3.1. Observations and Reduction 

We obtained deep images of NGC 7079 at the CTIO 0.9-m telescope on the night of 
1997 August 2-3 under photometric conditions. We used the TeK 2K#3 CCD in quad 
readout mode, with a gain of 3.3e~/ADU and read noise of 4.7e~ averaged over the four 
quadrants. The image scale in this configuration was (y.'396/pixel, for an unvignetted field of 
13'5 x 13f5. The exposure time for each image of the galaxy was 900 seconds; we acquired 
four exposures in /, two each in B, V, and R, and one exposure in U. The seeing through 
the night was ~ l'/5. For the /-band images, the telescope was moved between exposures 
so that the galaxy's image fell on different parts of the CCD, to allow the effects of bad 
pixels to be removed when the images were combined. The images were flat-fielded with 
twilight sky flats, and zero-subtracted as usual. Using the IRAF 7 task imalign, the J-band 
images were shifted to the nearest pixel, and then combined using the biweight, a robust 
statistical estimator of the mean that is insensitive to outliers (Beers et al. 1990). (Because 
of the small pixel scale relative to the seeing, integer pixel shifts were adequate, and avoided 
interpolation effects around image artifacts.) A constant sky value, obtained from image 
regions that were free of stars and galaxy light, was subtracted from each image. All were 
then cleaned of cosmic rays using IRAF's task cosmicrays. The resulting /-band image is 
shown in Fig. 1. A B — I map, also shown in Fig. 1, reveals little structure, apart from 
a nucleus bluer than the disk, indicating that there is no variable obscuration within the 
galaxy (as required for the TW method). 

Two standard stars were used for calibration, E9-47-U and E1-44-S (Graham 1982). 
We obtained the values of the Galactic extinction from Schlegel et al. (1998) via NED 8 . We 
did not measure the differential atmospheric extinction, but used the standard value from 
Hamuy et al. (1992). Since we observed NGC 7079 and the standard stars all at less than 
1.1 air masses, the resulting uncertainties are insignificant for our purposes. 



7 IRAF is distributed by NOAO, which is operated by AURA Inc., under contract with the National 
Science Foundation 

8 The NASA/IPAC Extragalactic Database (NED), is operated by the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under contract with the National Aeronautics and Space Administration. 
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3.2. Isophotal analysis 

After masking bad pixels, foreground stars, and all pixels outside 100" from the galaxy 
center, we fitted ellipses, with fixed centers, to single U, B, V and i?-band images, and 
to the combined /-band image, using the IRAF task ellipse. We determined integrated 
magnitudes using the fitted ellipses and IRAF's task bmodel. We do not correct these 
magnitudes for the galaxy's inclination. Our results are reported in Table 1, along with the 
total asymptotic U, B and V aperture-magnitudes given in RC3. Our magnitudes for NGC 
7079 are a bit fainter than RC3, possibly because we masked out the bright stars in the 
vicinity of the galaxy, and have not extrapolated the profiles to infinity. To measure the 
colors, we again used the task ellipse, but now constrained to measure the photometry 
in the same ellipses as were fitted to the /-band, so that magnitudes can be meaningfully 
compared in the same regions. The results are plotted in Fig. 2. Apart from the slightly 
bluer nucleus, the color profile is remarkably flat, with barely detectable color gradients. 

We fitted separate exponential disks to the photometry in B, V, R and / in the radial 
range 25" < R < 40", finding an average scale-length, Rd, of 17"1 ± 2"0 (corresponding 
to 3.1 ± 0.4 kpc at our assumed distance) and /io,s = 20.5 mag arcsec -2 . The latter value 
deprojects to 21.0 mag arcsec -2 , close to the Freeman (1970) limit (21.65±0.3 mag arcsec -2 ). 
At R ^ 45", the surface brightness profile steepens. This break radius remained unchanged 
when we fitted ellipses after changing the constant sky background level by ±la, and is 
identical for B, V, R and / bands, despite a difference in sky background level of three 
magnitudes. We conclude, therefore, that the break in the surface brightness is not an 
artifact of imperfect sky subtraction. Sharp disk truncations have been noted in photographic 
plate studies (van der Kruit 1979), while in a large sample of edge-on galaxies, Pohlen et al. 
(2002) found radial structure that is better fit by two exponentials, rather than by sharply 
truncated exponentials. His break radii (Rb/Rd = 2.5 ± 0.8) and surface brightnesses at Rb 
(fi R = 22.6 ± 0.6 mag arcsec -2 ) agree well with the values we obtain for the break in NGC 
7079 (2.6 and 22.0, respectively). While the fact that R b is identical in all bands argues 
forcefully that this is a real structure, it is clear that the disk profile continues past R b (cf. 
Narayan & Jog 2003). 

We use data at semi-major axis larger than 51"5 to calculate the median values of the 
disk ellipticity, e disk , and PA disk . Fig. 3 shows these results averaged over all colors. We 
find e disk = 0.36 ± 0.02 corresponding to an inclination i = 49?8 ± 1?7 for a razor-thin disk. 
If, however, the ratio of disk scale height to scale length is go — 0.2, we use the standard 
equation (Hubble 1926) cost = a/ (b 2 /a 2 — q$) /(l — g 2 ,), to obtain % = 51?3. We therefore 
adopt % = 50?6 ± 2?3, which is in good agreement with RC3 (i = 51?8). Our value of 
PA disk = 78?8 ± 0?1; the RC3 value for this parameter is 82°. 
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Since a bar represents an m = 2 perturbation of the disk density, a Fourier decomposition 
of the deprojected surface density is well-suited to measuring the semi-major axis of the bar, 
a B . Since two independent parameters, an amplitude and a phase, define each m ^ 
harmonic, two independent estimates of a B are possible from a Fourier decomposition. Our 
first estimate of a B , based on amplitudes, is given by the method of Aguerri et al. (2000), 
which has been used in other early- type barred galaxies (Debattista et al. 2002, ADC03). The 
method defines a B as the radius for which (I b /I ib ) = (I b /I ib )i ee 0.5[(I 6 / I ib ) max + (I b /I ib ) m ini 
where I b — 1 + A 2 , In, — 1 — A 2 and A 2 is the amplitude of the m = 2 Fourier component, 
normalized by that of m — 0. For NGC 7079 this prescription gives a B = 28'.'9. For a 
phase-based estimate of a B) we deprojected the ellipse fits and measured the phase angle of 
the resulting ellipses. The value of a B is then the largest radius out to which the phases are 
consistent with a constant, taking into account that deprojecting the bulge, which we did 
not subtract from our image for this analysis, results in a twist interior to the bar. We found 
a B = 21'.' 9 this way; we therefore adopt a B = 25'.'4 ± 3'.'5. Fig. 4 presents the bar-length 
analysis and also shows that Vw, the bar- angle in the intrinsic galaxy frame as measured 
from the major-axis of the disk, is ~ 58?5, which projects to PA bar = 32?8. 



4. Spectroscopy 

4.1. Observations and Reduction 

We observed NGC 7079 on 1997 July 28 with the CTIO 4-m telescope using the RFP 
imaging interferometer and the TeK 2K#6 CCD. The instrument has a circular field of view 
of radius ~ 85", and the detector was operated with 2x2 pixel binning, producing an image 
scale of / /70/pixel. The detector gain and read noise were 1.0 ±0.1 e~/ADU and 2.9 ± 0.2 
e _ , respectively. We used the "broad" RFP etalon, with a resolution FWHM ~ 4.2 A (150 
km s _1 ) and free spectral range of ~ 100 A at the Call line. The spectral instrumental profile 
of the RFP is well fitted by a Voigt profile, the convolution of a Lorentzian and a Gaussian 
of widths a i = 1.93 A and a g = 0.45 A, respectively. We used a blocking filter of FWHM 
100 A to select the desired spectral order. The wavelength calibration of the RFP was 
determined with a precision of ~ 0.04 A. Because the wavelength zero point and the optical 
center drift during the night (due to flexure, temperature changes, etc.), we took additional 
calibration measurements hourly through the night and incorporated these corrections in the 
data reductions. We found average drift rates of 0.086 A/hour and 0.36 pixels/hour for the 
zero point and the optical center, respectively. The temperature dropped during the night, 
until the telescope primary mirror was some 2.5°C warmer than the ambient air, producing 
poor seeing of l'/6 — 3'/0 (FWHM). There were also strong winds, reaching 50 mph by morning 
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and the night was non-photometric, with thin variable cloud. Nevertheless, we managed to 
obtain 24 exposures of 900 seconds each, spanning the wavelength range 8608 - 8631 A m 
steps of ~ lA (35 km s" 1 at 8618 A). We first stepped through the full wavelength range 
in steps of ~ 2A, and then filled in the intermediate wavelengths with a second scan (one 
wavelength sample was inadvertently repeated). 

The images were zero subtracted and flat-field corrected (with the median of 3 dome flats 
at each wavelength) in the the usual way. Cosmic rays were then cleaned by replacement 
using IRAF. In an imaging FP, the transmitted wavelength varies quadratically with the 
distance from the optical axis; for the RFP at Call the field is ~ 7.5 A bluer at the outer 
edge than on the optical axis. Thus the sky background in each frame varies radially due 
to the spectral energy distribution of the night sky. We offset NGC 7079 to the right half 
of the field of view, so that we could determine the sky spectrum from the left half of the 
field. Because some faint galaxy-light from the outer disk does reach to the left side of the 
frames, as well as a ghost image from reflections off the CCD and FP, we developed an 
iterative scheme to determine and subtract the sky. Briefly, the initial estimate of the sky 
spectrum was determined from the left half of the image, and subtracted from the entire 
image; this approximately sky-free image was used to estimate and remove the contributions 
of the outer galaxy and the ghost reflection from the left half of the original image, leading 
to an improved estimate of the sky spectrum. A second iteration of this procedure produced 
the final sky spectrum that was subtracted from the images. 

The images were then masked, shifted to a common frame and convolved with a radi- 
ally symmetric Gaussian to a common seeing. We did this twice, once for all 24 images, 
(seeing FWHM = 2'/6) and again for the best 18 images, excluding those obtained under 
poor conditions at the end of the night (resulting FWHM = 2'.'2). Since the atmospheric 
transparency varied by ~ 12% during the night, all images needed to be normalized to a 
common transparency. We determined the normalization coefficients by measuring the flux 
of seven bright stars in our field and assuming their mean spectrum over our wavelength 
range was flat and featureless. Finally, we generated noise maps which include photon, sky 
and read noise, and which account for the Gaussian and image registration smoothings. 

To determine the kinematic properties of NGC 7079, we fit a Voigt profile to the spectral 
data at each spatial position in the FP data cube. The data to be fit were smoothed with 
box car averaging using an adaptive kernel ranging from 1 x 1 to 7x 7 pixels based on the flux 
level, and the fits were weighted by the inverse square of the estimated flux uncertainties. The 
parameters of the fit were the continuum level, the line strength, the line-of-sight velocity, 
and the Gaussian velocity dispersion; maps of these quantities and their uncertainties were 
produced. We were able to obtain reliable fits within an ellipse of semi- major axis length of 
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~ 40" (pLi ~ 21 mag arcsec -2 ). We did not attempt to fit higher velocity moments, such 
as the Gauss- Hermite moments h 3 and h± (Gerhard 1993, van der Marel & Franx 1993). 
A comparison of the velocity fits produced from the 18 and 24 image data cubes showed 
no significant differences, mostly because of the poor quality of the last 6 images; in the 
following, therefore, we consider only the fits based on the better-seeing (18 image) data 
cube. 

Fig. 5 presents the resulting velocity and velocity dispersion maps. Other than a slight 
S distortion in the central regions due to the bar, the velocity distribution is that expected 
for an inclined, rotating disk. The velocity dispersion is largest in the nucleus and bulge, 
lower along the bar, and smallest in the disk. Examples of fitted profiles are presented in 
Fig. 6; these show successful fits from the galaxy center (panel 3) to fij = 21 mag arcsec -2 
(panels 1 and 7). By /// = 22 mag arcsec -2 , the signal-to- noise ratio is too low and only the 
noise is fit (panel 12). The spectrum of one of the foreground stars (panel 10) is essentially 
constant (with possibly a very weak Fe A8611 A line). As seen in panel 3, the FP scan did 
not extend far enough in wavelength to sample the continuum adequately in the high velocity 
dispersion regions in the center of the galaxy; thus while the velocities are well-determined 
here, the velocity dispersion measurements have larger systematic uncertainties. 



4.2. Comparison with slit spectroscopy 

Since this is the first application of the RFP to absorption-line spectroscopy of galaxies, 
it is important to confirm that our results agree with those obtained by traditional slit 
spectroscopy. For this purpose, E. M. Corsini kindly obtained and reduced for us a long- 
slit spectrum along the major axis of NGC 7079. The 1800 second exposure was taken 
in 16 November 2001 using the ESO New Technology Telescope (NTT) with EMMI in red 
medium-dispersion spectroscopic mode and using the grating No. 6 with 1200 grooves mm -1 
in first order with a l'.'O x 5'5 slit. The detector was the No. 36 Tektronix TK2048 EB CCD 
with 2048 x 2048 24/im pixels. It yielded a wavelength coverage between about 4840 A and 
5490 A with a reciprocal dispersion of 0.320 A pixel -1 . The instrumental resolution was 
1.19 A (FWHM) corresponding to a inst « 30 km s -1 at 5170 A. The spatial scale was 0'.'270 
pixel -1 . 

The spectrum was bias subtracted, flatfield corrected, cleaned of cosmic rays, corrected 
for bad pixels and columns, and wavelength calibrated using standard MIDAS 9 , as described in 
ADC03. The line-of-sight velocity, velocity dispersion, and the higher order moments, were 



9 MIDAS is developed and maintained by the European Southern Observatory 
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measured from the galaxy absorption features present in the wavelength range centered on 
the Mg line triplet (AA 5164, 5173, 5184 A), using the Fourier Correlation Quotient method 
(Bender 1990; Bender et al. 1994). 

In Fig. 7, we compare the slit spectrum velocity and dispersion data with data extracted 
from our RFP 2-D maps in a l'/O wide pseudo-slit along the major axis, smoothing along 
the pseudo-slit to match the slit spectrum smoothing. There is a zero-point offset of 14 
km s" 1 between the two velocity data sets, which has been added to the plotted data. The 
two data sets agree well, with very similar amplitudes and shapes of the rotation curves. A 
chi-square analysis suggests that the velocity uncertainties are under-estimated, and that an 
additional uncertainty of 8.4 km s" 1 should be added in quadrature to the error estimates. 
If all this additional uncertainty is assigned to the RFP data, the typical uncertainty of our 
velocity measurements is 12 km s -1 ; if the additional uncertainty is divided evenly between 
the two data sets, the typical RFP errors are 10 km s -1 (and typical slit spectrum errors 
are 8 km s _1 ). The RFP velocity dispersion data have more scatter and larger estimated 
errors than the slit velocity dispersions. As noted above, we expect the lack of sufficient 
wavelength coverage to affect the precision of the RFP dispersion measurements when the 
dispersions are large, as in most of the points along the pseudo-slit. Nonetheless, the two 
data sets are in good general agreement, tracing similar amplitudes and distributions of the 
velocity dispersion. Again, a chi-square analysis suggests that the uncertainties are under- 
estimated, and that an additional uncertainty of 20.1 km s _1 should be added in quadrature 
to the error estimates. If all this additional error is assigned to the RFP measurements, their 
typical precision is 31 km s -1 ; if the additional error is split equally between the two data sets, 
the RFP uncertainties are 27 km s _1 and the slit spectra errors are 15 km s _1 . There is some 
marginal indication that the velocity dispersions measured with the RFP are systematically 
larger than with the slit spectrograph. The mean difference of all the measurements is 8 
km s _1 , and the mean difference of the measurements at radial distances greater than 7" 
(where the RFP profiles are adequately sampled) is 15 km s" 1 ; these biases are smaller than 
the estimated uncertainties, and the spectral resolution of the RFP (corresponding to an 
instrumental dispersion of 64 km s -1 ) makes precision measurements of dispersions below 
50 km s" 1 problematic. 

Since the RFP and slit spectrograph data attain comparable results and accuracies, 
it is straightforward to compare the relative efficiencies of the two methods. Since the 
spectral resolution of this spectrograph was about 2.2 times better than that of the RFP, it 
should have been possible to widen the slit to 2" and still measure velocities with accuracies 
comparable to the RFP. In an 1800 second exposure, the spectrographic data reached to a 
limiting radius of 20", while the RFP data extended to 40", about 1.3 magnitudes fainter. 
The total RFP exposure time was 18 x 900 seconds. The sky and galaxy flux levels were 
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sufficiently high that both instruments were in the photon noise limit regime. Thus, taking 
into account the telescope apertures and the noise from the sky background, a single slit 
spectrogram with 2" slit would reach the same surface brightness about 3.3 x faster than the 
RFP. However, the RFP data extend out to ±20" along the galaxy minor axis; to cover this 
same area would require 20 spectrographic exposures. Thus the RFP is approximately 6x 
more efficient than a slit spectrograph for measuring the two-dimensional velocity field of 
this galaxy. 

We therefore conclude that the RFP data is in good agreement with the slit spectrum 
data, and that the RFP technique affords a reliable means for obtaining absorption-line 
galaxy kinematics with significantly greater efficiency that standard slit spectroscopy. 



Using the velocity field data, we obtained the stellar mean rotation velocity curve, V*(R), 
by fitting tilted-rings (Begeman 1987, Palunas 1996), with the galaxy center, inclination and 
PA disk fixed to the values found in the photometry. We started by determining the systemic 
velocity, V sys , obtaining V sys = 2670 ± 4 km s -1 . This is in good agreement with Lauberts 
& Valentijn (1989), who obtained V sys = 2670 ± 25 km s" 1 , and with the averaged values of 
Bettoni & Galletta (1997), V sys = 2678 ±35 km s _1 . This parameter was then also fixed and 
V*(R) determined using pixels at angle \cf>\ < 45° from the disk major axis in the intrinsic 
galaxy frame. The result is shown in Fig. 8. Subtracting V*(R) from the velocity map 
produced small residuals, mostly ^ 20 km s -1 except near the bright foreground stars. 

To determine the rotation curve of the galaxy we needed to correct for the asymmetric 
drift, a procedure which is somewhat uncertain. The asymmetric drift equation for the 
circular velocity, V c , is (e.g. Binney & Tremaine [1987] Eqn. 4-33) 



where p is the disk's volume density, and <j§ and &r are the tangential and radial velocity 
dispersions in the cylindrical coordinates of the galaxy's intrinsic plane. To proceed further, 
we needed to make a number of assumptions. Following ADC03, we assumed that the 
rotation curve is flat at large radii and that both <j b s and p decrease exponentially with 
radius. Thus the asymmetric drift equation becomes 



5. Kinematics and Bar Pattern Speed 



5.1. 



The rotation curve 






(3) 
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where a = <j z /<jr, <7 and R a are the parameters of the exponential fit to the observed major- 
axis velocity dispersion profile, and is the exponential scale-length of the surface density. 
As in ADC03, we use three values of a — 0.7, 0.85 and 1.0, appropriate for early-type galaxies 
(Gerssen et al. 1997, 2000). 

In Fig. 8, we present this asymmetric drift correction to points at R > 30". We used 
i?d = 17'.'1 from the photometry, as appropriate for the region of interest, and fitted an 
exponential to o" D b s along a 7-pixel- wide slit on the major axis. We found a circular velocity 
amplitude, V c fl at = 2561^ km s _1 . The radially- averaged mean streaming velocity in the 
same radial range is K = 200 ± 5 km s _1 . 

5.2. The pattern speed of NGC 7079 

We calculated the TW integrals for NGC 7079 in 9 strips in the ranges -14" <Y< -2'/3 
and 2" 3 < Y < 11".7. The limits on large \Y\ were chosen to avoid contamination by the 
bright stars flanking NGC 7079 near its the minor axis. The limits on small \Y\ were chosen 
to avoid a small (< 1 pixel) offset between the kinematic and photometric centers, to which 
V is sensitive because of the centrally peaked light distribution. In our standard case, we 
calculated the integrals inside an ellipse of semi-major axis, A max = 35", which does not quite 
reach the axisymmetric part of the disk (see Fig. 3), but where the surface brightness is more 
than five magnitudes below the central value (compare D03). We chose the function h(Y) 
for each integral to be the integrated surface brightness within the strip. We obtained error 
estimates for V from 100 Monte-Carlo realizations, where the observed velocity field was 
modified at each point by adding a velocity randomly drawn from a Gaussian distribution of 
width equal to the estimated velocity uncertainty at that point plus an additional 8.4 km s _1 
added in quadrature to account for the possible underestimate of the velocity uncertainty. 
In the standard case, we used for EpT, Y) the /-band map; to match the photometric and 
kinematic maps, we first convolved the former to the seeing of the latter, then used several 
of the stars in the field of both maps to compute a geometric transformation between them. 
This transformation required only a translation and a single scale for X and Y, with no 
rotation or shear. The results for the standard case are shown in Fig. 9; the best fitting 
slope is —6.5 ±0.1 km s _1 arcsec _1 (the negative sign here reflects only on our arbitrary 
choice of positive X). 

We checked that this value is not sensitive to various choices we made. We first tested 
that using the /-band photometry has not biased our result by using instead an image gen- 
erated by summing all the RFP frames. The resulting slope was —6.3 ± 0.3 km s _1 arcsec _1 , 
in good agreement with the standard case. We also checked that our limits on the integrals 
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are reasonable by confirming convergence of the values of X and V when we changed the 
semi-major axis, A max , of the ellipse within which the integrals were performed. In the range 
31'.'5 < A max < 38'.'5, the variations in V are all within the standard errors. We use the 
variations in X in 35" < A max < 38V5 to define the errors in X; these small errors lead to 
negligible change in the value of the fitted slope and its uncertainty. 

From this slope we calculate Q p = 8.4 ± 0.2 km s _1 arcsec _1 , and Dl = V c ^ t /Q p = 
30'.'6l{;3 (67%confidence interval). We then find that K = V cflat /(a B Q p ) = 1.2 ± 0.1 at the 
67% confidence interval and 1Z = 1.21q2 a t the 99% confidence interval, where we obtained 
the uncertainties by Monte Carlo experiments in which we varied V Ci fl at and ae uniformly in 
their error intervals and Q p assuming Gaussian errors. At a probability of higher than 94%, 
NGC 7079 has a fast bar. Note that, if the RFP over-estimates the velocity dispersions in 
the region of interest, then the real 1Z would be even lower; in the limiting (and unphysical) 
case that velocity dispersions are negligible in the region of interest, we obtain 1Z — 0.9 ±0.1. 



5.3. Variations in PA disk 

Our 2-D data also allowed us to study the effects of errors in PAdisk on the fitted slope. 
In the convention where the error in PAdisk, <W, moves PAdisk away from alignment with 
the bar when <5pa is positive (the same convention as in D03), we repeated our standard TW 
measurement as before but for —10° < 5 P a < 10°. The results, plotted in Fig. 10, reveal 
that the slope is quite sensitive to rather small 5pa- Errors in PAdisk of order 5°, for example, 
produce large errors in Q p , ~ 50% — 100%. Clearly, how accurately fl p can be measured 
depends strongly on how well PAdisk can be determined. Our high-precision determination 
of PAdisk (see Fig. 3) allows us to determine fl p with the small uncertainty quoted above. 
Note however that had we used the value of PAdisk from RC3, our value of Q p would have 
been some 35% larger. 

These results are consistent with those of D03 who studied AQ/Q p systematically by 
means of an A-body simulation. The sensitivity of Q p calculated here from the NGC 7079 
observations is, if anything, somewhat larger than in D03 (compare our Fig. 10 with Fig. 9 
in D03). Because PAdisk errors can also be produced by disk density ellipticities, D03 was 
able to show that such ellipticities had to be <^ 0.07 for the sample of SB0 galaxies with TW 
measurements. 
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6. Discussion and Conclusions 

We have shown that FP techniques can measure galaxy kinematics from stellar absorption- 
lines in early-type galaxies. To the best of our knowledge, this is the first successful demon- 
stration of this technique. Comparison with a traditional slit spectrogram along the major 
axis of the galaxy confirms the accuracy of the FP measurements. For the galaxy observed 
here, mapping the kinematics was about 6 times more efficient with the FP than with a slit 
spectrograph. For an object that fills more of the FP field of view, the advantage would be 
even greater. The fact that the spectral and spatial resolutions are decoupled with an imag- 
ing FP affords a great deal of flexibility in the data analysis, allowing extended measurements 
at lower spatial resolution in the faint outer regions of galaxies. Since the FP samples a large 
area on the sky, the determination and subtraction of the sky spectrum can be performed 
with high precision, again allowing measurement to fainter surface brightnesses. The full 2D 
format of FP kinematic maps permits unambiguous determination of the positions at which 
velocity samples are obtained, and allows optimal extraction of information in crowded fields, 
which can be difficult with slit data. 

There are some limitations to using a FP for absorption-line spectroscopy. The tempo- 
ral sampling of the spectrum limits the accuracy of reconstruction of the line profile; this is 
the reason why we did not pursue measurement of higher order moments of the velocity dis- 
tribution function. The next generation RFP, with a built-in transparency monitor, should 
help alleviate this difficulty The reduction of FP data has a reputation for difficulty; we 
believe that this is more due to the wealth of information available in a full 2D kinematic 
dataset than to the inherent difficulty of analyzing the FP data. Nonetheless, the techniques 
for handling FP data are not widely known, and FP software is neither as generally available 
nor as robust as the packages for slit spectroscopy. Finally, while we have shown the signif- 
icant efficiency of a FP for producing a full 2D map of the stellar kinematics of a galaxy, 
there are many classes of investigation for which such extensive data are not required. In 
such situations traditional slit spectroscopy is both simpler and more efficient. 

These problems should not, however, be overstated. The promise of FP absorption-line 
spectroscopy is demonstrated by the results obtained here. This is the tightest bound on 
1Z ever obtained on a barred galaxy using the TW method, with a 67% interval less than 
half that for the next best published case (ESO 139-G009, ADC03). More dramatically, the 
fractional uncertainty on fl p is less than 3%: typical uncertainties from slit spectroscopy are 
~ 20 - 30%. 

We showed that, in order to measure tt p to such high precision, the PA of the disk needs 
to be known quite accurately. Our's is the first demonstration on a real galaxy that the TW 
method is sensitive to small errors in PA. For NGC 7079, we were able to measure PA disk to 
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better than 1°. But if we had relied on the RC3 value of PA disk we would have obtained a 
value of tip ~ 35% larger. 

We found, at high confidence, a fast bar in NGC 7079. As in all previous measurements 
with good precision, the probability that 1Z — 1.7 is low: if we assume extreme values for a B 
and V^fl a t to maximize 71, we find that 1Z = 1.7 requires Q p = 7.3 km s~ 1 arcsec~ 1 , or over 5<r 
smaller than our measured value. Such fast bars continue to be a challenge for cuspy cold 
dark matter halos. 
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Fig. 1. — The J-band (left panel) and the B — I maps (right panel) of NGC 7079. North is 
left and west is up in each panel. The bright stars near the galaxy were used to normalize 
the RFP frames. The contours in the /-band map are spaced at 1 mag arcsec^ 2 , with the 
bold contour indicating \ij = 21 mag arcsec~ 2 . The only significant feature in the color map 
is the bulge at the center, with no evidence for a patchy obscuration. 
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Fig. 2. — The top panel shows the sky-plane surface brightness of NGC 7079 in the /-band. 
The sky level is indicated by the dotted horizontal line, while the dashed line indicates the 
best-fit exponential within 25" < R < 40". The radial profile exhibits a break at ~ 45" 
(fij ~ 21.2 mag arcsec" 2 ). The bottom panel shows the colors of NGC 7079. The squares 
are U — I, the circles B — I, the triangles V — I and the diamonds R — I. Note the quite flat 
colors (contamination by scattered light in the ?7-band is responsible for the feature seen in 
[/-/between 30" and 35"). 
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Fig. 3. — The sky-plane ellipticity (top) and PA (bottom) of NGC 7079. The data from the 
ellipse fits to the U, B, V, R and / bands have been averaged together for this figure. The 
disk's ellipticity and PA are determined from data with semi-major axis larger than 51". 5. 
This gives e d i S k — 0.36 ± 0.02 and PA disk — 78?8 ± 0?1. The average values of e disk and 
PA disk are shown by the dashed lines; the standard error on edisk is shown by the dotted 
lines, while the dotted lines in the bottom panel show PA disk ± 1°, illustrating that PA disk 
can be determined to better than 1° from these data. 
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Fig. 4. — The top panel shows the estimate of ob based on the amplitudes of the Fourier 
m = 2 deprojected moment. The dotted horizontal line shows the value of (h/Iib)i, from 
which as = 28'/9, indicated by the vertical dotted line, is obtained. The bottom panel 
presents the deprojected position angles relative to PA disk , ip, of the ellipse fits of Fig. 3 
(averaged over all bands), in the inner part of the galaxy. The vertical dotted line is at 
a B = 21'/9. The horizontal dotted line shows ^ bar ~ 58?5. 
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Fig. 5. — The line-of-sight velocity (left) and velocity dispersion (right) fields; north is left 
and west is up. In the velocity panel, contours are spaced by 50 km s -1 , centered on the 
systemic velocity and are solid (dashed) on the approaching (receding) side. The bold dashed 
lines in the velocity panel indicate the 2 regions within which we applied the TW analysis. 
Our limits on large \Y\ were chosen to avoid contamination by starlight, and the limits on 
small \Y\ are intended to avoid a small offset in the kinematic center. 
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Fig. 6. — Sample profile fits; each fit-panel (numbered in its upper-right and with the 
resulting V\ os indicated in its upper-left) points to the location in the galaxy (here shown 
by the /-band map from Fig. 1) for which it applies. The data points indicate the fluxes 
(in electrons) in the reduced data cube, while the green lines show our fitted profiles. The 
points are coded by the time order in which they were observed: the black (filled) points are 
the first 12 observed, while the blue (open) points are the last 6 observed. 
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Fig. 7. — A comparison of our data with slit data. The top panel shows line-of-sight veloc- 
ities, while the bottom one shows the velocity dispersions. The open symbols are slit data, 
the solid symbols are the RFP data. The error bars correspond to la deviations. 
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Fig. 8. — The observed stellar streaming velocities, V*, (top panel) and the asymmetric drift 
analysis for the amplitude of the rotation curve (bottom panel). The bottom panel shows 
the observed velocity dispersion along the major axis of the disk; the best fit exponential is 
indicated by the solid line. In the upper right corner of the lower panel are shown estimates 
of Vcflat (see Eqn. 2) for a = 0.7 (closed and open circles) at R > 30". The radial average 
for this a is indicated by the solid horizontal line. The dashed and dotted horizontal lines 
indicate the radial averages obtained by assuming a = 0.85 and a — 1.0, respectively. 
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Fig. 9. — The integrals X (top left) and V (bottom left) plotted as a function of Y . The bot- 
tom right panel plots V versus X\ fitting a straight line gives slope —6.5 ±0.1 km s _1 arcsec _1 
and x 2 — 2-2. Excluding the most discrepant point reduces x 2 to 0.7 without changing the 
slope or its error significantly. 
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Fig. 10. — The slope of the best fitting line (top panel) and the resulting error in JItw, 
the pattern speed derived from the linear regression of the TW integrals (bottom panel), as 
functions of error in disk PA, Spa- In the bottom panel, AQ/Q P = (fixw — ^p)/^V The 
solid circle in both panels represents the case when Spa = 0. 
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Table 1. Apparent magnitudes of NGC 7079. 



Filter 


m 




u 


13.2 


12.83 ± 0.13 


B 


12.6 


12.46 ± 0.13 


V 


11.7 


11.59 ± 0.13 


R 


11.1 




I 


10.4 





a RC3 (via NED) 

Note. — We estimate that our 
uncertainty is less than 0.1 mag- 
nitude in all bands. 



